Sphere packing with a geometric based compression algorithm

An efficient algorithm for the random packing of spheres can significantly save the cost of the preparation of an initial configuration often required in discrete element simulations. It is not trivial to generate such random packing at a large scale, particularly when spheres of various sizes and geometric domains of different shapes are present. Motivated by the idea of compression complemented by an efficient physical process to increase packing density, shaking, a new approach, termed compression algorithm, is proposed in this work to randomly fill any arbitrary polyhedral or cylindrical domains with spheres of various sizes. The algorithm features both simplicity and high efficiency. Tests show that it takes 181 s on a 1.4-GHz PC to complete the filling of a cylindrical domain with a total number of 26,787 spheres, achieving a packing density of 52.89%.     

Prediction of necking of high strength steel tubes during hydroforming—multi-axial loading

This article studies tubular hydroforming of high strength low alloy (HSLA) and dual phase (DP600) straight tubes under the action of end feeding loads. Experiments demonstrate that higher end feed loads enhance the formability of the tubes and increase the internal fluid pressure for onset of necking and bursting. Because of the action of the internal pressure and the axial compressive load, the onset of localization (necking) is due to a complex three-dimensional state of stress. Using free expansion experiments, approximate upper and lower bound strain-based forming limit curves are determined for the tube materials. These limit curves, in turn, are used to derive upper and lower bound extended stress-based forming limit curves [Simha et al., Prediction of necking in tubular hydroforming using an extended stress-based FLC. Transactions of the ASME Journal of Engineering Materials and Technology 2007;129(1): 36-47]. In conjunction with finite element computations that use solid elements to model the tube, these stress-based limit curves are used to predict upper and lower bound necking pressures under the action of end feed loading. These predictions of necking pressures, when an appropriate coefficient of tube-die friction is used, are found to bracket the experimentally measured necking pressures. Computations using plane stress shell elements to model the tubes are shown to give erroneous results, since the plane stress approximation is not valid when tubes are hydroformed in a die.     

磷矿中七种金属元素的光谱法测定

研究了磷矿样品中7种金属元素锰、钴、镍，锌、铬、钼及铋的分析方法，选择了灵敏度高干扰少的谱线，绘制出工作曲线，对样品进行了科学处理，不必进行化学分离即可对七种元素同时测定，经合成样分析、回收实验及样品测定，证明此方法精确可靠，其回收率在98％-102％之间。     

Bifurcation points and bifurcated branches by an asymptotic numerical method and Padé approximants

The aim of this work is to develop a reliable and fast algorithm to compute bifurcation points and bifurcated branches. It is based upon the asymptotic numerical method (ANM) and Padé approximants. The bifurcation point is detected by analysing the poles of Padé approximants or by evaluating, along the computed solution branch, a bifurcation indicator well adapted to ANM. Several examples are presented to assess the effectiveness of the proposed method, that emanate from buckling problems of thin elastic shells. Especially problems involving large rotations are discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

Second-order accurate constraint formulation for subdivision finite element simulation of thin shells

We present a new method for enforcing boundary conditions within subdivision finite element simulations of thin shells. The proposed framework is demonstrated to be second-order accurate with respect to increasing refinement in the displacement and energy norm for simply supported, clamped, free and symmetric boundary conditions. Second-order accuracy on the boundary is consistent with the accuracy of subdivision-based approaches for the interior of a body. Our proposed framework is applicable to both triangular and quadrilateral refinement schemes, and does not impose any topological requirements upon the underlying subdivision control mesh. Several examples from an obstacle course of benchmark problems are used to demonstrate the convergence of the scheme. Copyright © 2004 John Wiley & Sons, Ltd.     

A direct stiffness analysis of a composite beam with partial interaction

The use of the conventional semi-analytical stiffness method in finite element analysis, in which interpolation polynomials are used to develop the stiffness relationships, leads to problems of curvature locking when beam-type elements are developed for composite members with partial interaction between the materials of which it is comprised. The curvature locking phenomenon that occurs for composite steel–concrete members is quite well reported, and the general approach to minimizing the undesirable ramifications of curvature locking has been to use higher-order polynomials with increasing numbers of internal nodes. This paper presents an alternate formulation based on a direct stiffness approach rather than starting from pre-defined interpolation polynomials, and which does not possess the undesirable locking characteristics. The formulation is based on a more general approach for a bi-material composite flexural member, whose constituent materials are joined by elastic shear connection so as to provide partial interaction. The stiffness relationships are derived, and these are applied to a simply supported and a continuous steel–concrete composite beam to demonstrate the efficacy of the method, and in particular its ability to model accurately both very flexible and very stiff shear connection that causes difficulties when implemented in competitive semi-analytical algorithms. Copyright © 2004 John Wiley & Sons, Ltd.     

Reduced modified quadratures for quadratic membrane finite elements

Reduced integration is frequently used in evaluating the element stiffness matrix of quadratically interpolated finite elements. Typical examples are the serendipity (Q8) and Lagrangian (Q9) membrane finite elements, for which a reduced 2 × 2 Gauss–Legendre integration rule is frequently used, as opposed to full 3 × 3 Gauss–Legendre integration. This ‘softens’ these element, thereby increasing accuracy, albeit at the introduction of spurious zero energy modes on the element level. This is in general not considered problematic for the ‘hourglass’ mode common to Q8 and Q9 elements, since this spurious mode is non‐communicable. The remaining two zero energy modes occurring in the Q9 element are indeed communicable. However, in topology optimization for instance, conditions may arise where the non‐communicable spurious mode associated with the elements becomes activated. To effectively suppress these modes altogether in elements employing quadratic interpolation fields, two modified quadratures are employed herein. For the Q8 and Q9 membrane elements, the respective rules are a five and an eight point rule. As compared to fully integrated elements, the new rules enhance element accuracy due to the introduction of soft, higher‐order deformation modes. A number of standard test problems reveal that element accuracy remains comparable to that of the under‐integrated counterparts. Copyright © 2004 John Wiley & Sons, Ltd.     

Structural Appraisal of the Florentine Gothic Construction System

The nave of Santa Maria Novella, a Dominican church in Florence Italy is representative of a Florentine Gothic system of construction. This system, consisting of domical rib vaults on square nave bays, high side aisles, and crypto-buttressing, differs substantially from the French high Gothic system of even-level-crown rib vaults on rectangular bays, with flying buttresses over relatively low aisles. An investigation into the structural aspects of the Florentine Gothic construction system reveals that the domical vaults increase the longitudinal thrust and reduce the transverse thrust, which may require additional precautions to be taken during construction, but entails a less elaborate buttressing system. The domical vaults are found to perform well structurally, with an absence of tensile stresses, and are more suitable for a square bay. In general, the structural system of the nave of Santa Maria Novella is found to be the product of carefully considered structural design, which may be accepted as an alternative to the French Gothic construction system.     

A 3D contact smoothing method using Gregory patches

In this work, a method is developed for smoothing three‐dimensional contact surfaces. The method can be applied to both regular and irregular meshes. The algorithm employs Gregory patches to interpolate finite element nodes and provide tangent plane continuity between adjacent patches. The resulting surface interpolation is used to calculate gaps and contact forces, in a variationally consistent way, such that contact forces due to normal and frictional contact vary smoothly as slave nodes transition from one patch to the next. This eliminates the ‘chatter’ which typically occurs in a standard contact algorithm when a slave node is situated near a master facet edge. The elimination of this chatter provides a significant improvement in convergence behaviour, which is illustrated by a number of numerical examples. Furthermore, smoothed surfaces also provide a more accurate representation of the actual surface, such that resulting stresses and forces can be more accurately computed with coarse meshes in many problems. This fact is also demonstrated by the numerical examples. Published in 2002 by John Wiley & Sons, Ltd.